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Abstract - The record statistics of quantum standard map is shown to capture the classical 
transition to chaos. It is shown that in the mixed phase space regime the number of intensity 
records is a power law in the dimensionality of the state as opposed to the logarithmic growth for 
the random states. The exponent of this power-law is close to 0.5 at the critical value of the chaos 
parameter, K ~ 0.98 of the standard map, an exponent that is also obtained for random walks. 
These findings are based on the record statistics of complex, normalized random states for which 
we have shown that the probability of a record intensity is a Bernoulli process. 



Introduction. — The success of the many books and 
websites on record breaking events is a testimony of their 
enduring fascination. In an indexed data set, or discrete 
time series, a "record" is an event that is larger (upper 
record) or smaller (lower record) than all preceding ones 
[I] . In the following upper records are exclusively studied 
and are simply referred to as records. Questions about 
how the records increase with time, or the number of 
records set are of a natural interest in a variety of con- 
texts from sports to rainfall [2HS], and a mathematical 
theory of records for independent identically distributed 
(i.i.d.) random variables has been developed [7H9]. If 
{xt, t = 1, • • • N} is a finite time series, the first element, 
of the corresponding records series is x\ itself and at 
subsequent times t it will be Rt = max(xt, Rt-i). As Xt is 
a random variable, so is R t and properties of this random 
variable are of interest. 

The study of the statistics of extreme events plays a 
very important role in a variety of contexts from hydrol- 
ogy to climate change [T^. Apart from numerous statis- 
tical problems to which it naturally applies, the theory of 
extremes has also been used in the context of determin- 
istic dynamical systems [11]. One of the motivations of 
studying extreme values in quantum states is to under- 
stand whether deviations arising for dynamical systems 
originate from the localization of wavefunctions. This, in 
turn, is connected to the important issue of quantum er- 
godicity. The Tracy- Widom distribution [12] of the ex- 



treme eigenvalues of random matrices is also relevant to 
finding the fraction of entangled states [13]. In [14], the 
statistics of maximum and minimum intensities in com- 
plex random states, and, in the quantum standard map 
have been studied. However the study of records in ran- 
dom states and quantum dynamical systems has not yet 
been explored, and is the intention of this Letter. 

The statistical nature of quantum chaotic states have 
so far been explored mainly with the random matrix the- 
ory (RMT) [THOU], an d while indeed this has been suc- 
cessful, there remain many questions, especially for cases 
in which the classical dynamics is not fully chaotic, and 
also for many-body systems with few-body interactions. 
Deviations from universal RMT predictions may have in- 
teresting dynamical origins. Records is a relatively new 
statistical tool that is interesting because some of the fea- 
tures like number of records is "distribution free" for i.i.d. 
processes, that is are independent of the underlying prob- 
ability distribution for the variables, and therefore any 
deviations in this is a reflection of correlations. Large (or 
small) amplitudes is interesting from the point of view 
of localization, and therefore the study of extremes and 
records for wavefunctions is of interest. Indeed beyond 
quantum chaos, the theory of records can be of interest 
in the study of localization in general, including Anderson 
localization. 

Due to the complexity of the eigenfunctions of chaotic 
quantum systems, it is expected that the records in the in- 
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tensity would be distributed as those for random vectors. 
It may be noted that while the present work considers 
the intensities of the quantum states, this is interesting 
in itself as the statistical properties of the probabilities 
are associated with complete von Neumann position mea- 
surements. At the same time, the phases are not given 
a complete go by, as the momentum representation is si- 
multaneously studied. Indeed it is well-known that the 
momentum and position intensities together encapsulate 
essential quantum features such as the uncertainty princi- 

pleHU. 

Here general results on record statistics for a "time- 
series" employing the intensities of wavefunctions is ob- 
tained. These include results on interesting new quantities 
such as the number of records in a state, and the distri- 
bution of the position of the maximum. In the following, 
record statistics for normalized random states arc studied 
and subsequently the record statistics in the eigenstates 
of the quantum standard map - a paradigm of classical 
and quantum chaos - arc shown to reflect the nature of 
classical phase space, whether regular, mixed or chaotic. 
Significantly, not only is the transition to complete chaos 
reflected in such statistics, but the breaking of the last 
golden mean KAM tori is seen as a critical event in a way 
that will be elaborated. 

Record statistics for (5-correlated variables. 

Let \tp) be a normalized state vector, with compo- 
nents z n = (n\i/)) in a complete orthonormal basis \n). 
The "time-series" considered is simply x n = \(n\ip)\ 2 , 
n = 1, . . . , N, where N is the dimension of the Hilbert 
space. Thus the "time" in which the records are observed 
is not the real time, but could be some observable such as 
position or momentum. The question of records is then 
about statistics of large intensities, their positions, occur- 
ring along the index. Define the probability density for 
the record variable to be R, at time t as P(R,t). The 
average record (Rt) is given by the first moment. Let 
P(xi, . . . , xn) be the j.p.d.f. of TV random variables. The 
probability that the record at time t, Rt, is less than R is 
given by 



Q{R,t) 



rt 



dxi...dx t Pt(xi, ...x t ) 



(1) 



where Pt{x\, . . . , x t ) = J P{x\, . . . , XN)dxt+\-.-dxN is the 
marginal j.p.d.f. of the first t random variables. It follows 
that P(R, t) = dQ(R, t)/dR. 

Components of normalized complex random vectors 
Zl = (n\ip), have the j.p.d.f.: P{zi, z 2 , ■ ■ ■ , z N ) = 
(T(N)/ir N )6 fejLi \zj\ 2 - l) • This is also the distribu- 
tion of the components of the eigenvectors of the GUE 
or CUE (Gaussian or Circular unitary ensembles) random 
matrices. It is the invariant uniform distribution with re- 
spect to the Haar measure under arbitrary unitary trans- 
formations on the (2N — l)-dimensional sphere. The cor- 
relation induced by normalization becomes weak for large 



N. The intensities Xi = \zi\ 2 being the random variables, 
it is useful to define the j.p.d.f. as 



P(xi, ...,x n ;u) 



(2) 



where u is an auxiliary quantity - the actual j.p.d.f. cor- 
responding to u — 1. Defining 

Q(R,t;u)= / ~[dx k / P(x 1 ,...,x N ;u)dxt + i...dx N M 

Jo k=1 Jo 



leads to 



e- su Q(R,t;u)du 



T(N) 
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(4) 

Using the convolution theorem, and then setting u = 1 in 
Q(R, t; u) gives 

Q(R,t)= ^(-l) m QVl-mi?) Ar - 1 e(l-mi?), (5) 

m=0 ^ ' 

Hence ^P(i?,i) - ELiM) m+1 CX^ - - 
mR) N 2 0(1 — mR), the probability density that the 
record is R at time t. Note that the final record is also the 
global maximum of the sequence, hence the final-record 
statistics provides another way to calculate the distribu- 
tion of the maximum and it agrees with results obtained 
in [14] for random GUE vectors. A word of caution is 
however in order, only the last record and the global max- 
imum coincide with each other. It should be appreciated 
that the second record and the second largest member the 
sequence are not the same. Hence, record distribution is 
a different measure and should not be seen as a general- 
ization of the work in [T4] . 

The 5-correlated process studied here is a special 
case of the "product measure state" whose j.p.d.f. is 
P(x 1 ,x 2 , ...,x N ) = C N rj i= i f{xi)5(^ =1 Xi ~ M), where 
f[xi) is some function, M a mass parameter and Cn the 
overall normalization. These arise as the steady states 
in various dynamic models of non-equilibrium statistical 
physics such as the zero-range process for which record 
statistics is still unknown. However, it has been shown 
that for f(xi) ~ Ax" 1 , the one variable marginal dis- 
tribution undergoes a phase transition by increasing the 
density p (= M/N) for 7 > 2 [18]. This phase tran- 
sition from fluid phase to anomalous or normal conden- 
sate depending on the parameters 7 and p, reflects in ex- 
treme mass distribution too. The extreme value distribu- 
tion changes from Gumbel, via Frechet to a non-standard 
extreme value distribution as p changes from sub-critical 
via critical to super-critical values [19]. Thus the present 
study of records in a special case has wider motivations in 
such contexts as well. 

Returning to the problem at hand, it can be easily seen 
that the t th record is Gumbel distributed for large N, as 
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for large N and t 3> 1, Eq. ([5]) becomes 

Q(R,t) w (1 - cxp(-iVi?)) t w cxp(-tcxp(-iV J R)). (6) 

Since the Gumbel distribution is of the form 
exp[— exp[— (x — aAr)/&Ar]] where ajy and b]y are the 
shift and scaling, it follows that for the records statis- 
tics the relevant parameters are a at = \og(t)/N and 
bN = 1/N. The above form also appears in the limit 
when the correlations are ignored. 

The average record as a function of time is (R(t)) = 
1 — Jq Q{R, t) dR which is equal to 



N ^} ' m\m) N N fe' 



(T) 



where H t is a Harmonic number as defined above. Known 
asymptotics of the Harmonic numbers implies that 



(R(t)) = ^l 7 + Ht) 



1 oo 

2t ^ 



k=l 



B2k 

2kt 2k 



(8) 



where B^k are Bernoulli numbers, and 7 is the Eulcr- 
Maschcroni constant. Again, this presents a generaliza- 
tion of the average maximum intensity found in |14j which 
corresponds to t = N. 

Other interesting quantities include the average number 
of records, the distribution of the lifetimes of the records, 
the probability that the last record (which is also the max- 
imum) survives for a given time and so on. A remark- 
able well-known fact from the theory of records is that 
for i.i.d. variables these quantities are distribution-free, 
that is independent of the particular underlying distribu- 
tion p(x) pp. For example the average number of records 
(Nr) = Hm ~ log(iV) + 7 is indeed very small compared 
to the length N of the data set; typically records are rare 
events. These follow from a classic result [T][9] that the 
probability of a record occurring at position j is 1/j, inde- 
pendent of the past and future position of the records. In 
other words the probability of the position of the records 
is a Bernoulli process, Ber(l/j). 

As the rank order of elements distributed in an i.i.d. 
fashion docs not change due to a normalization induced 
delta function correlation, not surprisingly, the average 
number of records for random vectors (and also in case of 
product measure state) is, and can be proven to be, same 
as the (Nr) stated above. It is interesting to note that 
while the record values do depend, for finite N, on the 
global constraint, the number of records do not. Being 
distribution-free, the statistics of records is a direct mea- 
sure of correlations. There are few analytical results for 
correlated variables, one notable exception being a ran- 
dom walk where it has been shown that the number of 
records grows much faster as y/N while the probability 
that a record occurs at j decays as 1/ y/J [20] rather than 
the i.i.d. and random states results of log (AT) and 1/j re- 
spectively. In the case of the random walk the number 
of records is not a self-averaging quantity, the standard 
deviation being of the order of the mean. 



Records for eigenvectors of quantum standard 

map. — We turn from random vectors to eigenvectors 
of the quantum standard map. Classically, the stan- 
dard map on the unit torus is an area-preserving map: 
W ,p') = (q+P> P— (A'/27r) sin(27rg')) where as the param- 
eter K changes from zero to unity, the system undergoes a 
transition from integrablc (free rotor) to a mixed system. 
At about K sa 1, the last rotational KAM torus breaks 
allowing global diffusion in the momentum space |21| . If 
the standard map is unfolded to a cylinder, it displays 
normal diffusion in momentum for large enough K . When 
K 3> 5, the classical map is essentially fully chaotic. For 
such parameters this leads to quantum eigenstates which 
follow the CUE/GUE or COE/GOE results depending on 
the value of the phases a and /3, which control parity sym- 
metry and time-reversal symmetry respectively. One may 
think of j3 as originating from a threading magnetic flux 
line, while a originates from the potential itself. If (3 7^ 
and a 7^ 0, 1/2 we can expect that both the time-reversal 
symmetry and parity symmetry are broken and the typi- 
cal eigenstates would be like complex random states. For 
the record statistics, we calculate the statistical properties 
of eigenvectors of the Floquet operator [35] in the position 
basis of a kicked pendulum on a torus: U nn i = 



N-l 



A7^> xp 



x exp 



(m + /?) 2 . (m + fl , 
-in — h 2m — [n — n ) 



KN 
i— — cos 
2?r 



N 

2ix{n + a) 
N 



N 



(9) 



The dimensionality of the Hilbert space N is the inverse 
(scaled) Planck constant. Thus the "data" in this case are 
the various eigenfunctions and especially their intensities. 
Clearly this can, in general, depend on the space in which 
the eigenfunctions are represented. Thus for small values 
of K we expect there to be many localized states in the mo- 
mentum space while being nearly uniformly distributed in 
the position, and this will reflect in any studies of records 
or extremes. However for large K, position or momentum 
basis will be equivalent up to fluctuations. For all the 
numerical calculation, unless otherwise stated, N = 400, 
a = 0.25 and j3 = 0.25 are taken. 

The average record as a function of the index normal- 
ized by the dimension of the Hilbert space (which acts 
as "time" for these vectors) is plotted in Fig. [T] for var- 
ious values of K. While this agrees well in the chaotic 
region with the result on random states, there are signifi- 
cant deviations in the mixed phase space regime of K < 5; 
eventually, these disappear for K = 10. For K (= 0.3), 
in the position space, most of the records are set up by 
t/N = 0.5, originating in the very weakly broken par- 
ity symmetry. The momentum space average records are 
somewhat similar but mostly, they lie above the result for 
random states and are not affected much by the broken 
parity symmetry due to their localization. Indeed, the 
record for eigenfunctions of the quantum standard map 
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Fig. 1: The average record (R(t)), from the ensemble of eigen- 
states of the quantum standard map in the position represen- 
tation. The parameters used are N = 400 and K — 10 (highly 
chaotic), K = 5 (mostly chaotic), K = 1 (mixed phase space), 
and K = 0.3 (mostly regular). The analytical curve refers to 
the random state result in Eq. (f8j. In all cases a = /3 = 0.25. 



in the classically chaotic regime is found to be Gumbel- 
distributcd, see Fig. [2] Here, we have also plotted the 
distribution for the "record" when t = N, which refers 
to the maximum intensity, for which the explanation is 
found in [14]. For small N, deviations from the Gumbel 
distribution are evident. These are consistent with the 
exact result, P(R,t) derived above (shown in the inset of 

Fig. E]). 



- Exact 

-t = 300, N = 400 
-t = 400, N = 400 
-Gumbel 




Fig. 2: The distribution of the records when the index is t 
for eigenfunctions of the quantum standard map with K = 10. 
After rescaling and a shift, the distributions are of the Gumbel 
type, except for small N (see inset) where deviations are seen 
and the exact formula for P(R, t) is to be used. 

The distribution of the position of the maximum in the 
position representation is shown in Fig. [31 where one can 
see a transition to the uniform distribution along with the 
transition to classical chaos. The sharp peak at the center 
for small K (here K = 0.3) is interesting and deserves fur- 
ther comment. For small K, when there are many narrow 
classical resonances, a significant fraction of eigenfunctions 



that are localized on separatrices possess maximum inten- 
sity at or very close to q = 1/2. For instance, for N = 400 
when K — 0.1 and 0.3, about 75% and 50% of states are 
peaked at q = 1/2. As we increase K, more islands start 
appearing with turning points away from q = 1/2, and 
eigenstates localize on their interior, thereby the maxi- 
mum intensity shifts away from q = 1/2. This provides a 
qualitative mechanism leading to uniform distribution for 
maximum intensity for large K. The transition in clas- 
sical behavior from integrable (K = 0), to mixed phase 
space where islands of stability coexist within the stochas- 
tic sea (intermediate K), to fully-developed chaos is well 
captured by the quantity S]\[(m). In the inset of Fig. [3] is 
shown the maximum intensity of individual states (again 
A" = 400) vs their actual position. For K — 0.3, there are 
a large number of states which have maximum intensity 
at q = 1/2 (index = 200). This manifests in the form of 
the observed sharp peak. These get gradually destroyed 
on increasing AT; indeed such measures seem useful to pur- 
sue further in classifying states in the mixed phase-space 
regime. 
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*— K = [1, 1.5] 
*— K= [0.3, 0.8] 
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400 



Fig. 3: The distribution of the position of the last record set, 
which is also the maximum, for eigenfunctions of the standard 
map with N — 400 and for various values of K. The inset 
shows the maximum intensity /max of individual states and 
where it occurs. 

It has been shown above that for random A^-dimensional 
vectors there are on the average ~ log A" + 7 intensity 
records. While we can expect to see this for the standard 
map in the chaotic regime, the mixed and near-integrable 
regimes show a marked departure. The correlations lead to 
results that are similar to those for the random walks, with 
a power-law scaling in N. The ratio of standard deviation 
of the number of records to the average is found to be of 
order unity, indicating non-stationarity of the distribution. 
For very small K the number of records is simply O(N), as 
the eigenfunctions arc describable by smooth functions. In 
the mixed- regime, and, when K < 1, a power- law along 
with a logarithmic dependance is clearly indicated. In- 
triguingly, it is almost a pure power-law with exponent 0.5 
at the critical value K = 0.98. For moderate values such 
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as K = 2.3 a clear separation of the phase space into 
chaotic and regular regions and a subsequent separation 
of the quantum spectra seem to make the numerical fits 
unstable. Further, for very large K (say 9.8), the result 
from the random vector case applies (Fig. 2]). This re- 
flects the changing nature of the wavefunction intensities, 
and indicates that at criticality it is close to some kind 
oil/ f noise. Eigenvalue fluctuations in quantum chaotic 
systems are known to have f// 7 noise [23], and further 
work on the spectral properties of intensities is therefore 
called for. 
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Fig. 4: The average number of records in the eigenfunction 
of the quantum standard map of length TV vs N is plotted 
for various K values. Solid line (except for K = 9.8) are fitted 
expression a log(N)+bN s , for K = 9.8 the solid line is log(JV) + 
7. a, b, S for different K values are, (I) K = 0.3, a — 0.45±0.15, 
b = 0.306 ± 0.007, S = 1.001 ± 0.003,(11) K = 0.9, a « 0, 
b = 0.73 ± 0.06, S = 0.63 ± 0.01, (III) K = 0.98, a « 0, 
b = 1.2 ± 0.2, S = 0.50 ± 0.01 and (IV) K = 2.3, a = 0.9 ± 0.3, 
& = 0.7 ±0.6, 8 = 0.21 ±0.04. 

The average number of records as a function of K is 
presented in Fig. [5l where the effects of the scaling in the 
mixed regime manifest themselves. It is remarkable that 
while record numbers in both, the momentum and the po- 
sition representation, approximately converge after i(«2. 
it is only for K > 5 when there may be tiny islands (if at 
all) that it returns to the result of the random vectors. 
These observations establish the connection between clas- 
sical dynamics in all the regimes with the record statistics 
for the eigenvectors of the quantum standard map. Since 
the standard map is a paradigmatic model, we believe this 
connection to hold generally. 

Conclusions. — The probability that a record ap- 
pears in the intensities of correlated random vectors de- 
rived results on record of intensities of correlated random 
vectors at an index j is a Bernoulli process, the same as for 
i.i.d. variables. The quantum standard map presents itself 
as a generic instance where cigenfunctions become increas- 
ingly complex with the system parameter, K. Remark- 
ably, a signature of breaking of the last torus at K ~ 0.98 
is found to be in the exponent one-half in the power-law 
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K 

Fig. 5: An appropriately scaled and shifted average number 
of records {Nr) vs K for the eigenfunctions of the quantum 
standard map with TV = 400. 

obeyed by the number of records. It is very likely that the 
conclusions drawn on the quantum standard map would 
hold for other quasi-integrable and mixed systems. 

Various extensions are possible - for higher dimensions, 
the theory of records has a multivariate extension [TJ ; the 
role of simultaneous breakdown of parity and time-reversal 
on records of intensities would be worth studying too. 
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